Can a galaxy redshift survey measure dark energy clustering? 
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A wide-field galaxy redshift survey allows one to probe galaxy clustering at largest spatial scales, 
which carries invaluable information on horizon-scale physics complementarily to the cosmic mi- 
crowave background (CMB). Assuming the planned survey consisting of z ~ 1 and z ~ 3 surveys 
with areas of 2000 and 300 deg 2 , respectively, we study the prospects for probing dark energy clus- 
tering from the measured galaxy power spectrum, assuming the dynamical properties of dark energy 
are specified in terms of the equation of state and the effective sound speed c c in the context of an 
adiabatic cold dark matter dominated model. The dark energy clustering adds a power to the galaxy 
power spectrum amplitude at spatial scales greater than the sound horizon, and the enhancement 
is sensitive to redshift evolution of the net dark energy density, i.e. the equation of state. We find 
that the galaxy survey, when combined with CMB expected from the Planck satellite mission, can 
distinguish dark energy clustering from a smooth dark energy model such as the quintessence model 
(ce = 1), when c e ;$ 0.04 (0.02) in the case of the constant equation of state Wq — —0.9 (—0.95). An 
ultimate full-sky survey of z ~ 1 galaxies allows the detection when c c ^ 0.08 (0.04) for wo = 0.9 
(—0.95). These forecasts show a compatible power with an all-sky CMB and galaxy cross-correlation 
that probes the integrated Sachs- Wolfe effect. We also investigate a degeneracy between the dark en- 
ergy clustering and the non-relativistic neutrinos implied from the neutrino oscillation experiments, 
because the two effects both induce a scale-dependent modification in the galaxy power spectrum 
shape at largest spatial scales accessible from the galaxy survey. It is shown that a wider redshift 
coverage can efficiently separate the two effects by utilizing the different redshift dependences, where 
dark energy clustering is apparent only at low redshifts z ,$ 1. 
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I. INTRODUCTION 

Various cosmological probes such as supernovae in dis- 
tant galaxies Q, 0, the cosmic microwave background 
(CMB) sky 0, and the galaxy redshift surveys 0,H,IS0 
have given strong evidence that a dark energy compo- 
nent, such as the cosmological constant, constitutes ap- 
proximately 70% of the total energy density of the uni- 
verse, which derives the accelerating cosmic expansion 
at low redshifts. Because there is no plausible theoreti- 
cal explanation for its existence and magnitude (e.g. see 
0, 0), observational exploration of the nature of dark 
energy is one of the most important issues in modern 
cosmology as well as particle physics. 

An observational dark energy task we should first ex- 
plore to tackle this fundamental problem would be to 
determine whether the accelerating expansion is as a con- 
sequence of the cosmological constant. Relaxing this as- 
sumption leads to a generalized dark energy with dynam- 
ically evolving energy density, which can be characterized 
by a time-dependent equation of state w(a) = pae/Pde 
(the cosmological constant corresponds to w = — 1). The 
current level of accuracy in constraining the constant w- 
parameter is a{w) ~ 0.1 at ltr level with the best-fit value 
containing a model with w = — 1 (e.g. see 0, 0, Ho|). 
Future prospects aimed at pinning down the constraint 
on w by a factor of 10 have been extensively investi- 
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gated to address the usefulness of various cosmological 
experiments based on massive galaxy surve ys s uch as to- 
mographic weak lensing experiment (e.g. LlllT^|). the 
baryon oscillation experiment (e.g. [I4 fl4L Tl5|) and the 
cluster abundance experiment (e.g. [16]). 

Another important consequence of a generalized dark 
energy is the spatial perturbation, providing an indepen- 
dent clue to resolving the nature of dark energy from 
the equation of state. There are many previous efforts 
made to study how an inclusion of the dark energy per- 
turbations modify the cold dark matter (CDM) struc- 
ture formation scenarios: theoretical predictions on the 
modifications in the CMB power spectra and the galaxy 
power spectra [HlllllllliiiSillllllilill^and 
the observational exploration of the dark energy clus- 
tering signal from the WMAP data at low multipoles 
0- HE H3- US- The dark energy clustering is relevant 
only for structure formation at low redshifts z ,$ 1, where 
the net energy density is dominant to the cosmic expan- 
sion. In addition, a reasonable model including the dark 
energy perturbation predicts that the dark energy can 
cluster together with matter components at large spa- 
tial scales (inevitably at super-horizon scales), whilst the 
dark energy is smooth at small scales so that the addi- 
tional component does not largely change the small-scale 
structure formations such as galaxy formation. For these 
reasons, the dark energy clustering effect on the CMB 
observables is likely to appear only via the integrated 
Sachs- Wolfe effect (ISW) at low multipoles, where the 
Sachs- Wolfe effect generated at the recombination epoch 
is significant contamination to separate the ISW effect 
from the measured power spectrum. 
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The transition scale to divide the dark energy clus- 
tering and smooth regimes can be usefully modeled by 
the effective sound speed of dark energy j23, 113 • In 
this model, smoothness of dark energy can be tested by 
searching for the signature of the sound speed from cos- 
mological observables. Hu and Scranton [20] carefully in- 
vestigated a prospect of how the ISW effect measured via 
the angular cross-correlation between CMB and galaxy 
distribution can be used to probe the dark energy cluster- 
ing, assuming an all-sky, deep multi-color imaging galaxy 
survey out to z ~ 2. 

A galaxy redshift survey offers an alternative means 
for probing the dark energy clustering, through the mea- 
sured statistical properties of three-dimensional gravita- 
tional clustering at largest scales. Compared to the an- 
gular correlations, from a cosmological point of view, the 
redshift survey carries more information on the underly- 
ing mass distribution due to a gain of the modes along the 
line-of-sight or redshift direction. Therefore, the purpose 
of this paper is to, for the first time, investigate the ability 
of a galaxy redshift survey for testing the smoothness of 
dark energy from the measured galaxy power spectrum. 
In fact, there are several future plans for high-redshift 
galaxy surveys that are already being constructed or se- 
riously under consideration: the Fiber Multiple Object 
Spectrograph (FMOS) on Subaru telescope |3l[. its sig- 
nificantly expanded version, WFMOS |3J, the Hobby- 
Ebery Telescope Dark Energy eXperiment (HETDEX) 
|33| . and the Cosmic Inflation Probe (CIP) mission [34| . 
These surveys probe galaxies at higher redshifts z 0.5 
than the existing surveys such as SDSS and 2dF surveys. 
Such a high-redshift survey has several advantages over 
the lower redshift ones. First, given a fixed solid angle, 
the comoving volume in which we can observe galaxies 
is larger at higher redshifts than in the local universe, 
thereby reducing the sample variance. This would make 
it more straightforward to obtain a well-behaved survey 
geometry that can help measure largest-scale perturba- 
tions. Second, density perturbations at smaller spatial 
scales are still in the linear regime or only in the weakly 
non-linear regime at higher redshift, which gives us more 
leverages on measuring the shape of the linear power 
spectrum to break the parameter degeneracies. In this 
paper, we will consider the survey design close to the 
proposed WFMOS survey, which consists two types of 
surveys different in redshift coverage and survey area: 
0.5 < z < 1.3 with 2000 degree 2 and 2.5 < z < 3.5 
with 300 degree 2 , respectively. The redshift coverage of 
WFMOS is suitable to probe the dynamical dark energy 
whose effects are apparent only at low redshifts z £ 1 , as 
we will show below. 

The structure of this paper is as follows. In Sec. ITTlwe 
start with writing down background cosmological equa- 
tions, the Hubble expansion and the angular diameter 
distance, in terms of cosmological parameters. In Sec. lIIII 
we define the effective sound speed parameter to model 
dynamical properties of dark energy clustering, and re- 
view how the dark energy clustering leads to a scale- 



dependent modification in the linear power spectrum 
shape assuming the adiabatic initial condition. Sec. II VI 
defines the galaxy power spectrum in terms of the primor- 
dial power spectrum, the transfer functions and the scale- 
dependent growth rate of mass clustering. In Sec. we 
first define survey parameters intended to resemble a fu- 
ture survey being planned, and describe a methodology 
to model the galaxy power spectrum observed from a red- 
shift survey that includes the two-dimensional nature in 
the line-of-sight and transverse directions due to the cos- 
mological and redshift distortion effects. We then present 
the Fisher information matrix formalism that is used to 
estimate the projected uncertainties in the cosmological 
parameter determination provided the measured galaxy 
power spectrum. In Sec. I VII we show the prospects of the 
future survey for probing the dark energy clustering. In 
addition, we carefully study how a degeneracy between 
the dark energy clustering and massive neutrinos can be 
lifted by utilizing the redshift information of galaxy clus- 
tering. Finally, we present conclusion and some discus- 
sion in Sec. IV 111 



II. PRELIMINARIES: COSMOLOGY 

Throughout this paper, we work in the context of spa- 
tially flat CDM models for structure formation (e.g. see 
[35| and [36]). According to the Einstein general rel- 
ativity, the expansion history of the universe is given 
by the scale factor a(t), which is related to redshift 
via 1 + z = 1/a (we use a(to) = 1 today for our 
convention). The cosmic expansion during the matter 
dominated epoch is determined by density contributions 
from non-relativistic matter Sl m and dark energy i^de at 
present, in units of the critical density 3Hq/8ttG, where 
H = 100 h km s^Mpc" 1 is the Hubble constant. The 
Hubble expansion rate is given by 



H 2 (a) = H ( 



-3(l+iu ) 



(1) 



where we have assumed the constant equation of state of 
dark energy, 

Pdo(a) ldlnp do 

w = = -r-ji 1 = const. (2) 

Pde(a) 3 a in a 

Note that wq = — 1 corresponds to the cosmological con- 
stant. In this paper we restrict ourselves to a dark energy 
model with wq > — 1 for simplicity. 

The comoving angular diameter distance is expressed 
in terms of the Hubble parameter as 

. . r da 

DA{a) = J 1 mwW' () 

giving the distance-redshift relation via 1 + z = 1/a. 
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III. DARK ENERGY CLUSTERING FOR 
ADIABATIC INITIAL CONDITION 

To model dark energy clustering , w e employ a phe- 
nomenological model developed in [22]. In this model, 
the stress perturbation of dark energy, which governs 
properties of the dark energy clustering, is usefully speci- 
fied by the dark energy equation of state and the effective 
sound speed c e , where the latter is needed to model the 
non-adiabatic stress perturbation (also see ^3)- The ef- 
fective sound speed of a generalized dark energy is defined 



Spdc 



Spdc 



(4) 



in a "rest frame" coordinate system where the momen- 
tum density of the dark energy vanishes. For a more 
general coordinate system such as Newtonian gauge, the 
pressure perturbation of dark energy, which enters into 
the r.h.s. of the momentum conservation equation, can 
be expressed as 



5p dc = w Sp dc + Pdc(ce ~ w o) \$dc + 3 



a Udc 



2 c i o / 2 \ ac - 

= c e dp de + 3-(c e -Wo) — Pde, 
a k 



(5) 



where the first and second terms in the first line on 
the r.h.s. denote the adiabatic and non-adiabatic pres- 
sure perturbations, respectively, u dc denotes the pecu- 
liar velocity of dark energy, and we have assumed wo — 
constant in time. Note that setting Udc = in Eq. (J5J 
(corresponding to the rest frame of dark energy pertur- 
bations) reduces the pressure perturbation to the form 
given by Eq. (JIJ. From Eq. JSJ, one can find that pres- 
ence of the effective sound speed leads the pressure per- 
turbation to act as a restoring force against the gravita- 
tional instability just like the Jeans instability of baryon 
perturbation (e.g., see Sec. 16 in ^3)> m the overden- 
sity (<5d > 0) and underdensity (Sd c < 0) regions, the 
pressure perturbation prevents further collapse and ex- 
pansion, respectively, if the wavelength of perturbation is 
smaller than the sound horizon (see below). Conversely, 
for a model that has only the adiabatic pressure pertur- 
bation (5pd c = woSpdc), the pressure perturbation and 
the gravitational force in the momentum conservation 
equation become kwoSdc + (1 + w )k^ (^ is the gravita- 
tional potential), namely, the two terms have same sign 
because wq < and 'J cx — <5dc via the Poisson equation: 
therefore, the pressure perturbation accelerates the grav- 
itational instability on all scales. Note that, throughout 
this paper, we assume the constant sound speed with 
c c > and ignore the trace-free stress perturbation of 
dark energy for simplicity. 

We expect that this modeling to treat the sound speed 
as a free parameter can cover a broader range of dark 
energy models rather than working with a specific model 
of the Lagrangian of dark energy sector. For a scalar 



field dark energy model, the sound speed can be exactly 
computed in linear theory in terms of the kinetic energy 
as a function of the field 0, |2^. For the case of a 
canonical kinetic energy term such as the quintessence 
H [Hill E3, the sound speed c e = 1. For a more gen- 
eral modification of the kinetic term such as a phantom 
energy ,38) and a k-essence field [3^, ^(j j a general value 
of c c including c <C 1 is allowed. 

The dark energy sound speed sets a characteristic 
length scale in gravitational clustering: 



Ade,fs(a) 



da 



(6) 



/ " a*H{a) ■ 

We shall often call this scale the comoving sound hori- 
zon of dark energy at epoch a. One might define the 
corresponding wavenumber as fcdefs = 27r/Ad c ,fs- For the 
scalar-field dark energy model with canonical kinetic en- 
ergy (c c = 1), the sound horizon is comparable with the 
particle horizon scale. In this paper, we simply refer to 
a model with c e = 1 as a smooth dark energy model, 
because a galaxy survey we consider can probe only the 
galaxy clustering at scales inside the present-day hori- 
zon scale and there is no unique definition of smoothness 
of dark energy on super-horizon scales where the energy 
density varies depending on the freedom of time-slicing 
in general relativity. 

On large scales A > Ad e ,fsi the dark energy perturba- 
tion can grow together with non-relativistic matter com- 
ponents. In this case, the gravitational potential $ in 
Newtonian gauge 1 evolves in a flat universe (e.g., see 0) 
as 



$(fc,a) = 4> c {a)ni{k) 



where 



Ma) = 1 - 



H(a) 



(7) 



(8) 



and lZi(k) denotes the primordial curvature perturba- 
tion of wavenumber k. In the matter dominated regime, 
H(a) cx a~ 3 / 2 , leading to C — > 3/5. Here we have 
employed the adiabatic initial condition for dark energy 
perturbations; the initial conditions of all the compo- 
nents are set by the primordial curvature perturbation. 
It is worth noting that, for a cosmological constant model 
(wo = — 1), Eq. (JHJ can be rewritten as 



1 



H(a) 



da' 



1 



H(a) 



, da '2HZ(a>) 
H(a) r . , I 



H{al) 

a! dH 2 (a') 



da' 



o da '~HZ(a')a>3- 



(9) 



1 More precisely, <& is the curvature perturbation, but related to 
the potential perturbation via $ = — ^ when the trace-free stress 
perturbation is negligible, which is valid at redshifts of our in- 
terest. 
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where we have used the partial integral in the first line on 
the r.h.s., and dH 2 (a)/da = — 3H 2 n nM a~ 4 in the second 
line. The final expression is equivalent to the well-known 
formula for the growth rate of mass cluste ring , for exam- 
ple, given by Eq. (10.12) in (also see [Zif). 

On the other hand, on small spatial scales A < Ad e ,f s > 
the pressure perturbation of dark energy prevents clus- 
tering of the dark energy; the dark energy is thus consid- 
ered as a smooth component in this limit. The rcdshift 
dependence of the gravitational perturbation is given by 
$ oc g(a) via 1 + 2 = 1/q, and the growth rate g{a) can be 
computed (e.g. see [lj) by solving the following differen- 
tial equation with the initial conditions g(a m d) — 1 and 
dg I da(a m a) = in the deeply matter dominated regime 
a md (e.g., a md = 1CT 3 ): 



, d 2 g(a) 
din a 2 



+ [5 — 3ui(a)f2de(a)] 



dg(a) 



din a 

+3 [1 - w(a)] n de (a)g(a) = 0, (10) 



where the equation of state is generalized so that it is 
allowed to have a time dependence w(a), and fide(tt) is 
the dark energy density parameter at epoch a, defined 
as fide (a) = 87rGpde(a)/3i? 2 (a). For a cosmological con- 
stant model, the solution of Eq. (jl0|l can be given by the 
integral form of Eq. ©. Note that in the presence of 
the non-relativistic neutrinos the growth rate is further 
modified at scales below the neutrino free streaming scale 
(see 0), as we shall discuss below. 

Hence, introducing the effective sound speed of the 
dark energy enforces that the stress perturbation leads 
the dark energy perturbation to be gravitationally sta- 
ble on small scales so that small-scale structure forma- 
tions are not largely modified by the additional compo- 
nent compared to the concordance ACDM predictions. 



IV. SHAPE OF LINEAR GALAXY POWER 
SPECTRUM 



To model the linear power spectrum including the dark 
energy perturbation, we use the recipe in 22] to compute 
the transfer function of the potential perturbation ^>(k, z) 
assuming the adiabatic initial condition. Given the pri- 
mordial power spectrum, the power spectrum of the po- 
tential perturbation in the matter dominated regime is 
given by 



A%(k,z) = T 2 c {k,z)-5 2 n 



D cbv {k, z) 



T 2 {k) 



-l+n„+|a, ln(fc/fe ) 



,(H) 



where the primordial power spectrum, Pn(k) = (1Z 2 ) 7 is 
specified in terms of the primordial curvature perturba- 
tion S-ji, the spectral tilt n s and the running index a s , and 
Td c (fc, z) is the fitting formula given in 22] to describe the 
dark energy clustering contribution as explained below. 



Note that the primordial power spectrum shape is defined 
at fco = 0.05 Mpc" 1 . In Eq. dJ, we have imposed that 
the super-horizon potential perturbation is related to the 
primordial curvature perturbation as $(z m d) = (3/5)72.^ 
in the matter dominated regime. The growth rate for the 
total matter (CDM, baryon plus non-relativistic neutri- 
nos) perturbations are computed using the recipe devel- 
oped in [451 ] and normalized as -D c b„(fc,z) — > ag(z) at 
k — > 0. It is also noted that, throughout this paper, 
we employ the transfer function of matter perturbations, 
T(fc), with baryon oscillations smoothed out for simplicity 
(see Sec. IV. B in for the related discussion). 

According to the physical processes described in 
Sec. IIIII the fitting function Td e (k, z) in Eq. l(TT|) is given 
in terms of the growth rates at scales smaller and larger 
than the dark energy sound horizon as 



T de (k, z) 



1 + q 2 



3g(z)/[5Mz)]+q 2 '' 
where the variable q is defined as 



q — -Z—y Ade,fs(z)Ade,fs(2de), 
Z7T 



with the redshift Zd c being given by 



jOdc(^dc) _ j_ 
Pm(zdc) 71-' 



1 + Z dc 



fid 
fin 



3^0 



(12) 



(13) 



(14) 



Note that Td c — * 5<j) c {z)/[3g(z)] or 1 as k — > or oo, 
respectively, which reproduce the two asymptotic regimes 
discussed in Sec. lIIII For the cosmological constant model 
(wq = —1), Td c (k,z) = 1 because 3g(z)/[5<^ c ] = 1 as 
explained around Eqs. @ and l]10p. 

Galaxies are biased tracers of the underlying gravita- 
tional field. Hence, the linear power spectrum of galaxy 
distribution is related to the potential power spectrum 
via the Poisson equation as 



k 3 Pg{k lZ ) 

2tt 2 



= b\ 



( 2k 2 



\ 3H fi„ 



A|(fc,z) (15) 



where b\ is a scale- independent, linear bias parameter we 
assume. 

As carefully investigated in [23 (see Appendix A of 
that paper), the fitting formula (|12[) holds with batter 
than ~ 10% accuracy compared to the results obtained 
by directly solving the multi- fluid Boltzmann equations, 
for a range of the dark energy parameters we are in- 
terested in. In the presence of non-relativistic neutri- 
nos, however, it remains unclear how our treatment of 
Eq. I|ll|) holds accurate. Hence, an interesting issue is to 
develop an accurate transfer function for the total matter 
perturbations including the dark energy as well as neu- 
trino perturbations. This will be our future study and 
presented elsewhere. For our purpose, which is to esti- 
mate the ability of future surveys to probe dark energy 
clustering, our treatment of the galaxy power spectrum 
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FIG. 1: Upper panel: Shown is how the dark energy clustering 
amplifies the power of linear power spectrum P(k) at z = 0.6 
on spatial scales greater than the dark energy sound horizon, 
relative to the smooth dark energy model (P S m). Note that 
non-relativistic neutrinos are ignored in this plot. The dark 
energy equation of state is taken to wo — —0.8, and the sound 
speed is varied to show how the scale-dependent modification 
in P(k) changes with c c : with decreasing c c , the transition in 
P(k) amplitude appears at smaller length scales, i.e. larger fc. 
Shaded boxes around the curve of c e — 0.1 show the expected 
l-o" measurement errors for spherically averaged P(k) at each 
k bin for a survey covering 2000 deg 2 and 0.5 < z < 1.3 (see 
text for the details), but note that we have multiplied the error 
bars by a factor 1/2 for illustrative purpose. Comparing the 
solid curves with the error bars leads to a naive expectation 
that the dark energy clustering can be detected from this 
type of galaxy survey only when the sound speed is smaller 
than c e ~ 0.1. Middle panel: Redshift dependence of the 
power spectrum shape modification due to the dark energy 
clustering, where c c = 0.1 and too = —0.8 are fixed. Lower 
panel: As the fiducial value of wo is away from wo = — 1, the 
amplification in P(k) at small k is more enhanced, for a fixed 
Co = 0.1. In all the panels, the bold solid curve shows the 
same result for c e = 0.1, wo = —0.8 and z = 0.6. 



is accurate enough. This is partly justified by the fact 
that the following results are not largely changed with 
and without non-relativistic neutrinos, as will be explic- 
itly shown below. 

The top panel of Fig. illustrates how the dark en- 
ergy clustering induces a scale-dependent modification 
in the shape of linear power spectrum, P(k), at z = 0.6, 
compared to the smooth dark energy model. The dark 
energy equation of state is fixed to wq — —0.8, and 



the fiducial value of the sound speed c G is varied from 
c c = 1 to c = 10~ 2 . Note that c e = 1 roughly cor- 
responds to the smooth dark energy model such as the 
quintessence model, since the sound horizon is compara- 
ble with the Hubble horizon. It is clearly seen that, on 
spatial scales larger than the dark energy sound horizon 
given by Eq. J§}, the dark energy can cluster together 
with dark matter and thus enhances the power spectrum 
amplitude compared to the small-scale amplitude that 
matches the smooth model prediction (P/P sm = 1). The 
transition in P(k) amplitude appears from larger k with 
decreasing the sound speed c e , because the dark energy 
sound horizon gets shorter. It is worth noting that, since 
we have normalized the primordial power spectrum (see 
Sec. HV() . all the power spectra match at k <C fcdc.fs or 
k &de,fsj showing that the effect of dark energy per- 
turbation becomes independent of fc. Hence, one can 
measure c e only if the characteristic transition pattern in 
the power spectrum shape is accurately measured from 
observations. In other words, the power spectrum ampli- 
tude is not useful to constrain c e . 



The shaded boxes around the curve of c c = 0.1 repre- 
sent the 1-tJ measurement errors on P(k), reduced by a 
factor 2 for illustrative purpose, for the fiducial galaxy 
survey with redshift range of 0.5 < z < 1.3 and survey 
area of 2000 deg 2 corresponding to the comoving volume 
of 3.6/i -3 Gpc 3 . Note that the errors are for the power 
spectrum spherically averaged over angle. It is appar- 
ent that the future survey of our interest allows a precise 
measurement of the linear power spectrum shape, achiev- 
ing a few % accuracies in each fc-bin at fc <L 0.05h Mpc -1 , 
while the errors are dominated by the sample variance at 
smaller fc. A quick look at this plot leads to a naive 
expectation that, only if the sound speed is sufficiently 
small such as c c £ 0.1, one could measure the transition 
pattern in P(k) from the galaxy survey. In the following, 
we shall carefully study how well a future galaxy survey 
can probe the dark energy clustering fully taking into 
account degeneracies between cosmological parameters. 



In the middle panel we show that the dark energy per- 
turbation more amplifies the power spectrum amplitude 
at the large spatial scales as one goes to a lower redshift, 
because the net dark energy density is more prominent at 
lower redshifts. Observationally this means that a galaxy 
survey of lower redshift z ^ 1 is more suitable for exploit- 
ing the dark energy perturbations than high-z ones. 



Finally, the lower panel of Fig. U shows how the effect 
of dark energy perturbation varies with the fiducial value 
of the equation of state wq. With increasing wq from 
wo = — 1, the dark energy perturbation more enhances 
the amplitude of P(k) at the large scales, because the 
dark energy density keeps prominent at higher redshift. 
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0.17 
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0.81 
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1.0 


0.19 
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1.55 
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0.21 


2000 


1.2 


0.5 


1.7 


Gz3 (2.5 < 


z < 3.5) 


3.0 


0.53 


300 


1.2 


1.0 


3.3 



TABLE I: Galaxy survey specifications that we assume in this paper. We consider two types of galaxy surveys, named as 
"Gzl" and "Gz3", intended to resemble the WFMOS fiducial survey |32|. The former sees galaxies at z ~ 1 with a fixed sky 
coverage of 2000 deg 2 , while the latter probes z ~ 3 galaxies with 300 deg 2 . V B and n g are the comoving survey volume and 
the comoving number density of sampled galaxies for each redshift slice, respectively. Note that V a is computed for our fiducial 
cosmological model with Q m — 0.27, Qdc = 0.73 and Wo = —0.8. z cen ter denotes the centering redshift of each redshift slice, 
and fc max is the maximum wavenumber below which information in the linear power spectrum can be extracted. (We do not 
use any information above fc max in the Fisher information matrix analysis.) "Bias" denotes the assumed linear bias parameters 
of sampled galaxies. 



V. METHODOLOGY 



A. Survey Parameters 



To derive a realistic parameter forecast, we employ the 
galaxy survey parameters that are chosen to resemble the 
future surveys that are under serious consideration. As 
shown below (see Eq. |18|L the statistical error of the 
galaxy power spectrum measurement is limited by the 
survey volume, V s , as well as the mean number density 
of galaxies, n g . There are two advantages for a higher 
redshift survey over the current surveys probing the uni- 
verse at z ^ 0.3. First, given a fixed solid angle, the co- 
moving volume in which we can observe galaxies is larger 
at higher redshifts than in the local universe, thereby re- 
ducing the sample variance. In addition, it would be rel- 
atively straightforward to obtain a well-behaved survey 
geometry, e.g., a cubic geometry that would be helpful to 
probe a largest-scale galaxy clustering as well as handle 
the systematics under control. Second, density fluctua- 
tions at smaller spatial scales arc still in the linear regime 
or only in the weakly non-linear regime at higher redshift, 
which gives us more leverages on measuring the shape of 
the linear power spectrum. 

We employ the survey parameters that match the fidu- 
cial survey design of WFMOS .32] , consisting of two types 
of redshift surveys different in redshift coverage and sur- 
vey area: 

• Gzl: 0.5 < z < 1.3 and D, s = 2000 deg 2 

• Gz3: 2.5 < z < 3.5 and D, s = 300 deg 2 

where the names Gzl and Gz3 stand for the "Ground" - 
based galaxy surveys probing the universe at z ~ 1 and 
z ~ 3, respectively 

Because we have limited knowledge of how galaxies 
have formed within the CDM hierarchical clustering sce- 
nario, it is of critical importance to figure out an optimal 
survey strategy of which type of galaxies (emission lines) 



are targeted to achieve the desired scientific goals, given 
the spectrograph specifications (sensitivity, the number 
of fibers, wavelength coverage, e tc.) and the number of 
nights allocated. It was shown in 01 (also see |46|l that a 
survey having n g P g ^ 3 over the range of wavenumbers 
considered is close to an optimal design from a cosmo- 
logical, not practical, point of view, and the hypothetical 
WFMOS survey parameters were defined in this regard. 
In this pap er, we simply adopt the survey parameters 
used in [13|, which are summarized in Table [I] 



B. A Galaxy Power Spectrum in Redshift Space 

We employ the method developed in p| (also see 0) 
to model a galaxy power spectrum observable from a red- 
shift survey. The power spectrum measures how clus- 
tering strength of galaxy distribution varies as a func- 
tion of 3-dimensional wavenumber, k (or the inverse of 
3-dimensional length scale). However, we cannot know a 
true position of a galaxy in real space from the observ- 
ables, angular position and redshift. Rather, we have 
to assume a fiducial cosmology to convert the observed 
angular position and redshift of galaxies into positions 
in 3-dimensional space. Since this fiducial cosmology 
generally differs from the true cosmology, we could in- 
troduce distortion in the inferred distribution of galax- 
ies. This cosmological distortion effect is the so-called 
Alcock-Paczynski (AP) effect Taking into account 

the effects of the cosmological distortion, linear redshift 
distortion and galaxy bias, the galaxy power spec- 
trum in the linear regime can be expressed in terms of the 
real-space galaxy power spectrum P g (k,z) (see Eq. |15j ) 
as 



P s (k 



'fid_L, Kfid| 



) 



D A (z)l d H(z) 
D A (z) 2 H(z) Rd 

xP g (k,z), 



1 + 13 



(16) 
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with 



k± = ^ } & , d fcfidJ 

Da[z) 



h = 



H(z) 
H Rd {z) 



Midll 



(17) 



where k = (k\ + fc 2 ) 1 / 2 , fey and k± denote the com- 
ponents of the wavenumber parallel and perpendicular 
to the line-of-sight direction, respectively, and (3 denotes 
the linear redshift space distortion taking into account 
the non-relativistic neutrino effect (see [44j), defined as 
(3 = -(l/bi)dlnD chu (k,z)/dln(l + z). D A (z) and H{z) 
are the comoving angular diameter distance and Hub- 
ble parameter, respectively, and the quantities with sub- 
script 'fid' denote the quantities in the fiducial cosmo- 
logical model. Note that P g (k,z) depends on the scale- 
independent, linear bias parameter b\ as P g oc b\ (see 
Eq. P3). 

The redshift and cosmological distortion effects are ex- 
tremely powerful for constraining cosmological parame- 
ters. Because the structure formation scenario predicts 
that the power spectrum P(k) has characteristic features 
such as the broad peak from the matter-radiation equal- 
ity, scale-dependent suppression of power due to baryons 
and non-relativistic neutrinos, the tilt and running of the 
primordial power spectrum, and the baryonic acoustic os- 
cillations, the distortion effects can be precisely measured 
with future galaxy surveys from the distorted features 
and help break degeneracies between parameter determi- 
nation quite efficiently. There are notably two examples 
for the prospects. First, the baryon oscillation peaks in 
P{k) can be used as a standard ruler, allowing precise 
measurements of H(z) and Da(z) to constrain the equa- 
tion of state of dark energy jl3l IT3. ITU EtSj . Second, 
measuring the redshift distortion helps break strong de- 
generacy between the power spectrum amplitude and the 
galaxy bias, allowing a precise determination of the neu- 
trino mass Ull. 



C. Fisher Matrix Analysis 

In order to investigate how well one can constrain the 
cosmological parameters for a given redshift survey, one 
needs to specify measurement uncertainties of the galaxy 
power spectrum. When non-linearity is weak, it is rea- 
sonable to assume that observed density perturbations 
obey Gaussian statistics. In this case, there are two 
sources of statistical errors on a power spectrum measure- 
ment: the sampling variance (due to the limited number 
of independent wavenumbers sampled from a finite sur- 
vey volume) and the shot noise (due to the imperfect 
sampling of fluctuations by the finite number of galaxies). 
To be more specific, the statistical error on measurement 
of P{k) at a given wavenumber bin ki is given in [49| by 



AP s (ki) 



Ps(ki 



2(2tt) 



V;fc 2 AfcA^ 



f 



f 



n qPs {ki 



(18) 



where n g is the mean number density of galaxies, V s is 
the comoving survey volume, and /i is the cosine of the 



angle between k and the line-of-sight direction. Note 
that we have assumed that the galaxy selection function 
is uniform over the redshift slice we consider and ignored 
any boundary effects of survey geometry for simplicity. 

The first term in the bracket on the r.h.s of Eq. i|18|) 
represents sampling variance. Errors become indepen- 
dent of the number density of galaxies when sampling 
variance dominates (i.e., P s 3> n g over the range of k 
considered), and thus the only way to reduce the errors 
is to survey a larger volume. On the other hand, the 
second term represents shot noise, which comes from dis- 
creteness of galaxy samples. When shot noise dominates 
(P s <C fi g ), the most effective way to reduce noise is 
to increase the number density of galaxies by increas- 
ing exposure time per field. Note that for a fixed n g 
the relative importance of shot noise contribution can be 
suppressed by using galaxies with larger bias parameters, 
bi, as P s oc b\. 

We use the Fisher information matrix formalism to 
convert the errors on P s (k) into error estimates of model 
parameters |l3i l50j ] . The Fisher matrix is computed from 



Faf3 



dfi 

1 Jk 

n g P s (k,/j) 



k 2 dk 



dlnP s (k,/i) d\nP s {k,n) 



dp a 



n g P s (k,fi) + 1 



(19) 



where p a expresses a set of parameters. The partial 
derivative with respect to a given parameter p a is eval- 
uated around the fiducial model. When combined with 
the CMB constraints, we simply add the CMB and galaxy 
fisher matrices, — F^ v + F™ B , to obtain the joint 
constraints. The la error on p a marginalized over the 
other parameters is given by a 2 (p a ) = (i r,_1 ) aQ , where 
F^ 1 is the inverse of the Fisher matrix. We follow the 
method described in 12] when we consider projected con- 
straints in a two-parameter subspace to see how the two 
parameters are correlated. 

To calculate F a p in Eq. ifH?|) . we need to specify k min 
and fc max for a given galaxy survey. We use the up- 
per limit, fcmax, to exclude information in the non- linear 
regime, where the linear theory prediction of density fluc- 
tuations, Eq. (IS) , becomes invalid. Following we 
adopt the values of fc max for each redshift slices as given 
m Table [Q As for the minimum wavenumber, we use 
kmin = 1CP 5 Mpc -1 , which gives well-converged results 
for all the cases we consider. 



D. Model Parameters of a CDM Model 

The Fisher matrix formalism assesses how well given 
observables can distinguish the true ("fiducial") cosmo- 
logical model from other models. The parameter fore- 
casts we obtain depend on the fiducial model and are 
also sensitive to the choice of free parameters. We include 
a fairly broad range of the CDM dominated cosmology: 
the density parameters are f2 m (= 0.27), il m h 2 (= 0.14), 
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and S!b/i 2 (= 0.024) (note that we assume a flat universe); 
the primordial power spectrum parameters are the spec- 
tral tilt, n s {= 1), the running index, a s (= 0), and the 
normalization parameter of primordial curvature pertur- 
bation, Siz(= 5.07 x 10~ 5 ) (the values in the parentheses 
denote the fiducial model). The linear bias parameters, 
b\ , are included for each redshift slice as given in TablelD 
The dark energy parameters, the equation of param- 
eter wq and the sound speed c e , are allowed to vary in 
order to study how the constraint on c e changes with 
the assumed fiducial values. We also include the non- 
relativistic neutrinos as free parameters (see ^ij for the 
details). The neutrino oscillation experiments have pro- 
vided strong evidence that two or three of standard three- 
flavor neutrinos have finite mass (e.g. see |H3 for a 
review). We assume throughout this paper that some 
of the three species are massive and have become non- 
relativistic until the present. The non-relativistic neu- 
trinos affect the structure formation by suppressing the 
growth of matter density fluctuations at small spatial 
scales owing to their large velocity dispersion. Hence, the 
non-relativistic neutrinos induce a scale-dependent mod- 
ification in the power spectrum shape. For a range of the 
total neutrino mass implied from the neutrino oscillation 
experiments, the effects of non-relativistic neutrinos and 
dark energy clustering could be degenerate in a galaxy 
power spectrum, especially when the dark energy sound 
horizon is comparable with the neutrino free-streaming 
scale. We shall below study this issue quantitatively. The 
neutrino parameters we assume are the ratio of the total 
neutrino mass density, f v {= = 0.01), correspond- 

ing to the total neutrino mass 0.13 eV, and the number of 
non-relativistic neutrino species, N^(— 2). The fiducial 
values of j v and iV" r are consistent with the neutrino os- 
cillation experiment results (see Fig. 3 in 0]) as well as 
the upper bounds set from the cosmological experiments 
0, 0, IE3 • In summary, for a survey which consists of N s 
redshift slices, we have 10 + N s parameters in total: 

p a = {Qm,Qmh ,&bh , 8n,n s , a s , w , c e , f v ,N™ 

M*i)>'--M*tt.)}- ( 2 °) 

As shown in the literature (e.g., see Q), a galaxy 
survey alone cannot determine all the cosmological pa- 
rameters simultaneously, but would leave some param- 
eter combinations degenerated. Combining the galaxy 
survey constraints with the constraints from CMB tem- 
perature and polarization anisotropies can be a power- 
ful way to lift parameter degeneracies (5^. When com- 
puting the Fisher matrix of CMB, we employ 7 parame- 
ters: 6 parameters (the parameters above minus the dark 
energy, neutrino and bias parameters) plus the Thom- 
son scattering optical depth to the last scattering sur- 
face, t(= 0.10). Note that our treatment is rather con- 
servative: we do not use any CMB constraints on the 
dark energy parameters and the non-relativistic neutri- 
nos. The Fisher matrix for the joint experiment is given 
by adding the CMB Fisher matrix to the galaxy Fisher 
matrix as F a[3 = F^ + F^ B . We entirely ignore the 
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FIG. 2: A significance of discriminating the dark energy clus- 
tering model from a smooth dark energy model (c e = 1) as a 
function of the fiducial value of c c , expected from the galaxy 
survey combined with Planck. The fiducial value of the dark 
energy equation of state is kept fixed to wq — —0.8. The bold 
solid curve shows the forecast for the nominal survey design 
of WFMOS (see Table Q, while the dot-dashed curve shows 
the result when the Gzl survey alone is used. The dashed 
curve shows the result for an ultimate all-sky survey of z ~ 1 
galaxies. 



contribution to the CMB from the primordial gravita- 
tional waves. We use the publicly-available CMBFAST 
code |54j to compute the angular power spectra of tem- 
perature anisotropy, C; TT , i?-mode polarization, C EE , 
and their cross correlation, C ; TE . Specifically we con- 
sider the noise per pixel and the angular resolution of 
the Planck experiment that were assumed in [HjJ. Note 
that we use the CMB information in the range of mul- 
tipole 10 < I < 2000, and therefore we do not include 
the ISW effect at low multipoles Z 10 which the dark 
energy perturbation may affect. 



VI. RESULTS 

A. Forecasts for Testing Dark Energy Clustering 

As can be seen from Fig. ^ t ne effect of dark energy 
clustering on the linear power spectrum is only apparent 
for variations that span an order of magnitude in c c (also 
see [3(j for the similar discussion). We therefore use 

Pc c = logioCo, (21) 

as a free parameter in the Fisher matrix calculation 
rather than directly using c e . Then, we define the fol- 
lowing quantity to estimate the significance of how well 
the future galaxy survey we consider can discriminate be- 
tween a given c e model and a smooth dark energy model 
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with c P = 1: 



where cr(log 10 c e ) 



logi 



cr(log 10 C c ) 



(22) 



(F- 



the 68% C.L. error on 



310 - Y V- /pc c Pc ( 

p Co marginalized over the other parameter uncertainties. 

The bold solid curve in Fig. [2 shows the forecasted sig- 
nificance (S/N) 2 to discriminate between a dark energy 
clustering model and the smooth model, for the WFMOS 
fiducial survey (Gzl plus Gz3 in Table H} combined with 
Planck. Note that the fiducial value of wq = —0.8 is 
assumed here. It is clearly seen that the dark energy 
clustering can be distinguished at significance (S/N) > 1 
as the fiducial value of c e gets smaller than c e ~ 0.08, 
where the dark energy sound horizon scale enters into 
a window of the wavenumbers where the galaxy power 
spectrum can be precisely measured with the galaxy sur- 
vey (see Fig. QJ. However, for too low sound speeds 
such as c ^ 0.01, S/N saturates at c c 0.01 (or the 
dashed and dot-dashed curves even have a peak around 
c c ~ 0.01). This is because the transitions in P g (k,z) 
occur at k comparable with or larger than the maximum 
wavenumber /c max above which the linear power spectrum 
information is not available due to the non-linearities in 
structure formation. On the other hand, the significance 
drops off rapidly as c — > 1 because the transitions oc- 
cur at too large spatial scales where the galaxy power 
spectrum measurement is largely limited by the sample 
variance. One may also notice that the curves appear 
to have a jagged feature around c = 0.05. This arises 
from degeneracies between dark energy clustering and 
non-relativistic neutrinos. This can be partly seen from 
the thin solid curve, which shows the results when the 
neutrino parameters f v and -/V" r are fixed to the fidu- 
cial values; the jagged feature disappears. We shall come 
back to this issue below in more detail. 

The dot-dashed curve shows the result when the Gzl 
survey alone is combined with Planck. Comparing the 
solid and dot-dashed curves clarifies that the constraint 
on dark energy clustering is mainly from the lower red- 
shift survey compared to the z ~ 3 survey, because the 
net dark energy density is more prominent at lower red- 
shifts (see Fig. The dashed curve shows the forecast 
for an ultimate full-sky redshift survey of z ~ 1, hypo- 
thetically obtained by enlarging the survey area of z ~ 1 
by a factor of 20. In this case, the model with c c ^ 0.2 
can be detected at more than a 1-a level. This result 
might be compared with Fig. 8 in [3(j; the authors inves- 
tigated the ISW effect to probe dark energy clustering 
via the CMB-galaxy cross-correlation measurements as- 
suming a full-sky CMB map and a full-sky, multi-color 
galaxy imaging survey out to z ~ 2. Interestingly, our 
results show a similar-level precision of constraining dark 
energy clustering to theirs, although our method does not 
include the CMB information on dark energy clustering 
(see Sec. lVDf) . and the survey area (/sky ~ 0.05) and the 
redshift range are much smaller and narrower than theirs, 



"ft 



X 

a 
S 



0.2 



0.1 



0.05 



0.02 



0.01 










— • 




^ — 














- / 




/ 




/ 




/ 




/ ^ 




/ 


■ 






/ / 




1 / 




1 / 
/ / 

■ ■ ■ ■ ' 


— Gzl + Gz3 





-1 



-0.95 



-0.9 



-0.85 



-0.8 



-0.75 



'O.fid 



FIG. 3: A maximum dark energy sound speed with 
(S/N)\ Co > 1 against the fiducial value of wo. The solid and 
dashed curves are the results for WFMOS (Gzl plus Gz3) 
and the full-sky z ~ 1 survey, respectively. As the underly- 
ing true cosmology approaches to the cosmological constant 
model, one can detect dark energy clustering from a galaxy 
survey only when the sound speed c e is sufficiently small. 



respectively. This would be partly because a galaxy 
power spectrum contains more independent modes of lin- 
ear fluctuations in the line-of-sight direction compared to 
the angular correlations. It is also noticed that compar- 
ing the dashed curve with the dot-dashed (or solid) curve 
shows about factor 20 increase in (S/N) 2 , corresponding 
to the gain in the sky coverage / s ky This implies that the 
constraint on c is mainly from the galaxy survey, not as 
a result of breaking the parameter degeneracies obtained 
by combining the galaxy survey with CMB, as expected. 
In summary, we conclude from the results in Fig. [21 that 
a wide-field galaxy redshift survey by its own offers an 
alternative, vital opportunity for testing the smoothness 
of dark energy. 

Obviously, the ability to detect dark energy clustering 
from a galaxy survey depends upon the net dark energy 
density or equivalently upon the equation of state wq in 
our setting for a fixed c e , provided f2dc and f2 m are tightly 
constrained. Fig. [3] shows a maximum sound speed de- 
tectable with (S/N) > 1 as a function of the given fiducial 
value of wo ■ As dark energy becomes close to the cosmo- 
logical constant (wq = —1), it is getting more difficult 
to detect dark energy clustering: dark energy clustering 
can be discriminated only if the sound speed is sufficiently 
small so that the scale-dependent transitions in P g (k) are 
measurable to within the measurement accuracies. To be 
more specific, the model with wq ^ — 0.96 allows a test 
of the dark energy smoothness for the WFMOS survey, 
while the model with wq £ — 0.98 can be tested for the 
ultimate full-sky survey. Alternatively, if the true cosmol- 
ogy has wo = —0.9 (—0.95), the dark energy clustering 
can be detected when c e 0.02 (0.04), while the full-sky 
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survey allows the detection when c c ;$ 0.04 (0.08), a fac- 
tor of 2 improvement over WFMOS. 



B. Parameter Degeneracies 

To more quantitatively understand the forecast for 
constraining the dark energy clustering shown in the pre- 
ceding subsection, it is useful to examine how the param- 
eters are degenerate with each other given the survey 
data sets. Table Hi) gives the correlation coefficients be- 
tween the dark energy sound speed and the other param- 
eters for the WFMOS survey (Gzl plus Gz3) combined 
with Planck. In the context of the Fisher information 
matrix formalism, the correlation coefficient is defined as 



Pc c Pc 



(23) 



where p a denotes one of the model parameters. When the 
coefficient |r| = 1, the two parameters are totally degen- 
erate, while r = means they are uncorrelated. While 
the coefficients varies with the fiducial model, we have 
assumed w = —0.8 and c c = 0.05 for the fiducial values 
of the dark energy parameters, where c G = 0.05 roughly 
corresponds to the dip scale of the curves in Fig. [2] (also 
see below for a more detailed discussion). Table Ull shows 
that the dark energy sound speed is most degenerate with 
the non-relativistic neutrino parameter iV" r as expected, 
because a change in N™ affects the shape of linear power 
spectrum via the change in the neutrino free-streaming 
scale 2 , but is insensitive to the P(k) amplitude; the sup- 
pression rate in P{k) amplitude due to the neutrino free- 
streaming is primarily sensitive to another neutrino pa- 
rameter /„ as AP(k]/P(k) « — 8/„ in the small-scale 
limit (see Fig. 2 in |44()- On the other hand, the sound 
speed is only weakly correlated with the other param- 
eters. Encouragingly, therefore, while the parameters, 
O m /i 2 , ftbh 2 , n s and a s , also modify the shape of power 
spectrum, combining the galaxy survey with Planck is a 
very powerful wa y to efficiently lift the parameter degen- 
eracies (also see |44| for the extensive discussion on the 
degeneracies between the neutrino parameters and the 
primordial spectrum shape parameters). 

It is also worth noting that the sound speed is only 
weakly correlated with the parameters controlling the 
galaxy power spectrum amplitude such as 6n and b±, 
again verifying that the constraint on c e is primarily from 
the shape of galaxy poser spectrum rather than the am- 
plitude information. Our method given in Sec. IV Bl in- 
cludes the redshift distortion effect on the galaxy power 
spectrum, which significantly helps break the degeneracy 



2 Note that the neutrino free-streaming wavenumber is propor- 
tional to the absolute mass scale of non-relativistic neutrino 
species and therefore depends o n th e neutrino parameters of our 
interest as k u>ts oc } V /N™ (see Q). 
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FIG. 4: Shown is how each redshift slice of the WFMOS sur- 
vey contributes to the dark energy constraint shown in Fig. [5] 
The number in the parenthesis for the label of each curve de- 
notes the comoving survey volume for each redshift slice (see 
Table|IJ. The lowest redshift slice (0.5 < z < 0.7) yields most 
contribution even though the slice has the smallest survey vol- 
ume and the smallest fc max , as dark energy is more prominent 
at lower redshifts. Also interesting is a jagged feature appears 
around c c « 0.05 when one redshift slice alone is used, while 
the jagged feature is significantly weakened when all the slices 
are combined. The reasons for this to happen are explained 
in text as well as the subsequent figures. 



between the power spectrum amplitude and the galaxy 
bias. We have checked that ignoring the redshift distor- 
tion effect little changes our results. Hence, constraining 
c e with the galaxy survey, if detected, is generally robust 
against uncertainties involved in galaxy bias. 



C. Dark Energy Clustering and Massive Neutrinos 

In this subsection, we investigate in more detail the 
degeneracies between the dark energy sound speed and 
the non-relativistic neutrinos. 

In Fig. 0] we start with looking into the significance 
(S/N) 2 , as in Fig. El but for the case when only one of the 
redshift slices listed in Table [I] is used. All the curves for 
each slice of Gzl show a more prominent jagged feature 
around c c = 0.05 than the solid curve for all the slices 
combined, implying that the constraint on c e is weakened 
due to the parameter degeneracies, because the dark en- 
ergy sound horizon becomes comparable with the neu- 
trino free-streaming scale for a redshift range accessible 
from the galaxy survey (we shall show this more explic- 
itly below). However, most important result shown here 
is the parameter degeneracies can be quite efficiently bro- 
ken, when all the redshift slices are combined, as shown 
by the bold solid curve. This implies that a survey with 
wider redshift coverage can efficiently separate the two 
effects from the measured P(k). 
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TABLE II: The correlation coefficients between the dark energy sound speed parameter p Cc (see Eq. 1211 1 and the other parameter 
p a for the parameter error estimates for WFMOS combined with Planck. The fiducial values for dark energy parameters are 
set to wq — —0.8 and c e = 0.05. The sound speed is most degenerated with the non-relativistic neutrino parameter iV" r . 
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FIG. 5: Effects of the dark energy clustering (solid curves) 
and the massive neutrinos (dashed) on the linear power spec- 
trum shape, relative to the smooth dark energy model. The 
five curves from top to bottom are the results for z — 
0.6, 0.8, 1.0, 1.2 and 3, respectively. We assume c c = 0.05, cor- 
responding to the jagged feature in the curve of Fig. [I] and 
wo — —0.8 for the dark energy fiducial model and assume 
f v — 0.01 and iV" r = 2 for the neutrino mass and the num- 
ber of non-relativistic neutrino species. Note that the dashed 
curves for the neutrinos are multiplied by an arbitrary factor 
so that P v (fc)/Pdc.sm — » 1 at large k, taking into account the 
fact that the galaxy power spectrum amplitude involves un- 
certainty related to galaxy bias. It is clear that the redshift 
dependences of the two effects are quite different. The two 
arrows in the horizontal axis denote fc max for the z ~ 1 and 
z ~ 3 surveys as given in Table [I] 



Figs. and [5] give a more quantitative explanation on 
the results in Fig.Q] Figure |S] compares the effects of the 
dark energy clustering and the non-relativistic neutrinos 
on the linear power spectrum shape, for 5 redshift epochs. 
We have assumed wo = —0.8 and c c = 0.05 for the fidu- 
cial values of dark energy parameters as in Table [HI It 
is intriguing to find the scale-dependent modifications in 
P g (k,z) due to these two effects are quite different: the 
dark energy clustering induces sharper fc-transitions in 
P g {k, z) as well as stronger redshift dependence than the 
non-relativistic neutrinos do. This difference can be ex- 
plained as follows. The neutrino free-streaming effect 
on P g (k, z) has accumulated over a longer time duration 
since the neutrino became non-relativistic, at z m w 240 
for m v ,i » 0.13 eV of our fiducial model, resulting 
in the modification of the linear power spectrum shape 
over a wider range of k. One the other hand, the dark 
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FIG. 6: Projected 68% ellipses in f v and c e sub-space. The 
error ellipses shrink more as more redshift slices are added. 
In particular, for the nominal WFMOS survey, a z ~ 3 sur- 
vey helps break the degeneracy between the two parameters 
efficiently, even though the slice is insensitive to dark energy. 
Note that the prior a(wo) = 0.05 is employed. 

energy clustering affects P g (k,z) only at low redshifts 
z £ 1. One can thus utilize the different redshift- and 
fc-dependences to separate the two effect from the galaxy 
power spectrum, if wider redshift and wavenumber cov- 
erages are available. 

Fig. H3 shows the projected 68% C.L. contours in the 
f u -c e sub-space, obtained when only one of the redshift 
slices for WFMOS is used or some/all of the slices are 
combined. Note that the prior a(wo) = 0.05 was em- 
ployed in this plot 3 . It is clear that, as a wider redshift 
coverage is available, the degeneracy between /„ and c e 
determination can be more broken. It is interesting to 
note that adding the z ~ 3 slice into the z ~ 1 slices can 
further shrink the error ellipses, even though the dark 
energy clustering is negligible at z ~ 3. This is because a 
higher redshift survey is more powerful for constraining 
the neutrino parameters due to a gain in the linear fluctu- 



3 We have so far ignored the baryon acoustic peaks in P(k), and in- 
cluding the peaks is powerful to constrain the dark energy equa- 
tion of state. To take into account this, we assume the prior 
0"(uiq) = 0.05 taken from Il3l . 
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FIG. 7: Difference in erg when a generalized dark energy con- 
tribution is included and the power spectrum is normalized by 
the primordial curvature perturbation. When wq > — 1 and 
the sound speed is large enough, erg is reduced compared to 
the ACDM prediction, because the structure formation slows 
down due to the more accelerated cosmic expansion at inter- 
mediate redshifts. On the other hand, when the sound speed 
is sufficiently small, the dark clustering enhances the power 
spectrum amplitude, yielding greater erg. The plot shows 
c e ~ 0.005 is the critical value of c c , where the two effects 
above cancel out so that the resulting erg is independent of c e 
and wo, and similar to the ACDM prediction. 

ations down to larger k (i.e. larger fc max ), allowing a more 
precise measurement of the suppression in the small-scale 
power of P{k) due to the neutrino free-streaming. 



spectra for the CMB-based methods, as we have so far 
employed (see Sec llVfl . Then, a s can be computed once 
the linear power spectrum shape and the growth rate 
from z ~ 1000 to present are specified to within the 
measurement accuracies (e.g., see |57| for a fitting for- 
mula of OS for a flat CDM model with dark energy of 
constant equation of state) . For this procedure, assum- 
ing a model with dynamical dark energy could affect an 
estimate of as compared to the cosmological constant 
model because the generalized dark energy could either 
suppress the growth rate due to the accelerated expan- 
sion or add the dark energy perturbation to the power 
of total mass perturbation. Which of these opposite two 
effects is dominant to determine a% depends on whether 
the dark energy sound horizon scale is greater or smaller 
than 8ft. -1 Mpc. Hence, if a discrepancy between the 
a 8 values inferred from the CMB and the local universe 
observations is found, it may be as a result of the dark en- 
ergy contribution, and it would be worth to keep having 
attention to explore such a signature from the cosmolog- 
ical probes at different redshifts. 

Fig.[7|shows how the value of ag changes as a function 
of c e and Wo, where other cosmological parameters are 
kept fixed to the fiducial values. When the sound speed 
c c is large enough compared to 8/i _1 Mpc, a s is reduced 
with increasing wq from wq = — 1 by slowing down the 
growth rate of mass clustering. On the other hand, if 
the sound speed becomes sufficiently small, albeit unre- 
alistic, the dark energy clustering increases as as a result 
of adding a power to the mass clustering at scales rele- 
vant for erg . Another interesting result is that the critical 
sound speed c GlC r ~ 0.005 emerges where the resulting as 
remains unchanged compared to the ACDM prediction. 



D. Effect of Dark Energy Clustering on erg 

The cosmological parameter as, the rms mass fluctua- 
tion today in spheres of radius 8/i _1 Mpc, is one of the 
most important parameters for characterizing the clus- 
tering strength in large-scale structures at low redshifts 
(e.g., |55|). However, for the moment cr 8 is relatively less 
accurately constrained compared to the other parame- 
ters within the concordance CDM model. The three-year 
WMAP data alone determines a$ to about 7% accuracy 
Q. In addition, methods based on independent data 
sets such as CMB, weak lensing and Ly-cv data sets show 
slightly inconsistent best-fit values of ag, for example, as 
implied in Fig. 7 in |3| (also see [56|): the WMAP favors 
slightly lower as than the methods looking at the local 
universe. 

Since the CMB observables are primarily sensitive to 
the perturbations at z ~ 1000, the present-day fluctua- 
tions, as, are rather an output parameter derivable from 
the parameters constrained by the CMB spectra. It is 
becoming common to use the primordial curvature per- 
turbation amplitude 5-jz to normalize the linear power 



VII. CONCLUSION AND DISCUSSION 

In this paper we have investigated the ability of a fu- 
ture galaxy redshift survey to test the smoothness of dy- 
namical dark energy, which is another important conse- 
quence of a generalized dark energy with w(z) ^ — 1 and 
provides independent information on the nature of dark 
energy from that carried by the equation of state. The 
dark energy clustering signatures can be measured via 
a scale-dependent transition in the power of the galaxy 
power spectrum appearing at scales comparable with the 
dark energy sound horizon (see Fig. It was shown 
that, for WFMOS survey, the sound speed can be de- 
tected at more than a 1-a level, if the sound speed is small 
enough as c e < 0.04 (0.02) when w = -0.9 (-0.95) (see 
Figs. |21 and 01. An effective way to improve the ability is 
to enlarge the survey volume especially for low redshift 
slices at z ~ 1. An ultimate full-sky survey could im- 
prove the lower bound on the detectable c e by a factor of 
2. 

Another interesting possibility of the future galaxy sur- 
vey is the use of the galaxy power spectrum to weigh the 
neutrino mass, as investigated in |44j |. We carefully in- 
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vestigate possible degeneracies between the dark energy 
clustering and the non-relativistic neutrinos in the galaxy 
power spectrum (see Fig. and also see for the re- 
lated discussion on the degeneracy between the neutrino 
mass and the dark energy equation of state). We showed 
that having a wider redshift coverage can efficiently sep- 
arate the two effects by utilizing the different redshift de- 
pendences; the dark energy is prominent only at low red- 
shifts z ^ 1. In addition, modified gravity theories have 
received much recent attention as an alternative expla- 
nation of cosmic acceleration without dark energy, where 
gravity is weaker than the Einste in g ravity at scales com- 
parable to the horizon (e.g., |^ |6l|, |63 ) . An impor- 
tant test to discriminate this possibility from the dark 
energy model is to simultaneously explore the cosmic ex- 
pansion history and the growth of structure formation 
from cosmological experiments, because the two possibil- 
ities predict different growth rate even for an identical 
expansion history. Once again, wider redshift coverage 
of a galaxy redshift survey would be powerful to make 
a robust test of discriminating between effects of dark 
energy clustering, massive neutrinos, and modified grav- 
ity, which warrants a detailed study (also see [f^] for the 
discussion on the interplay between the massive neutri- 
nos and the modified gravity in the shape of mass power 
spectrum) . We hope that the results shown in this paper 
will give a useful guidance to designing a future survey if 
the dark energy clustering is desired to pursue as one of 
the science goals. 

We have assumed throughout this paper that dynami- 
cal properties of the dark energy are specified in terms of 
the equation of state and the effective sound speed, with 
the adiabatic initial condition. A merit of this approach 
is we need not assume a specific form of the Lagrangian 
of dark energy sector. However, the limitation is there is 
no guarantee this modeling can be applied to a full range 
of generic dark energy models. For example, we have 
ignored a contribution from the trace-free stress pertur- 
bation of dark energy, which is valid if the dark energy 
is a scalar field (e.g. see |2j| for the extension). Since we 
have little idea of what dark energy is, it is worth explor- 
ing a more general modeling of dark energy properties 
and investigate the resulting effect on structure forma- 
tion. Unless such a theoretical model is available, we 
cannot extract all the cosmological information inherent 
in the future high-precision cosmological data sets in an 
unbiased way. 

Another non-trivial assumption we have made in this 
paper is the adiabatic initial conditions for dark energy 
perturbations, without any rigorous reasoning. If the 
dark energy field had existed since the inflationary epoch 
in the early universe, the dark energy perturbations nat- 
urally contained iso-curvature modes, and the effect just 
emerges at low redshifts. Hence, a galaxy redshift survey 
could by its own open up a new window for exploring 
the iso-curvature modes in structure formation to probe 



the physics in the early universe, complementarity to the 
CMB or the CMB-galaxy cross correlation [HES]. Fur " 
thermore, if the primordial non-Gaussianity (e.g. see |6^| 
for a thorough review) is dominantly imprinted onto the 
sector of dark energy, the galaxy survey allows us to ex- 
plore the signal via, e.g., the bispectrum measurement of 
galaxy clustering. These are our future study and will be 
presented elsewhere. 

The linear perturbation theory makes secure predic- 
tions on structure formation, which allows an accurate 
interpretation of the cosmological data sets from the de- 
tailed comparison. A main obstacle that a galaxy red- 
shift survey could contain in this procedure is uncertain- 
ties involved in the galaxy biasing, where unknown, non- 
linear gastrophysics may pollute the linear information 
even at large scales. Although we have assumed a scale- 
independent, linear bias, we do not think that a more re- 
alistic bias such as a scale-dependent bias is a significant 
contamination to the dark energy clustering constraint. 
Because the dark energy clustering induces characteristic 
scale-dependent modification and strong redshift evolu- 
tion for P(k) (see Fig. [I}, these properties can be very 
robust to discriminate the dark energy clustering signals 
from the galaxy bias contamination or more generally 
other systcmatics. 

Finally, we comment on the non-linear gravitational 
clustering, which could also provide a contamination to 
the linear theory predictions. Encouragingly, however, 
the non-linear effects can be to some extent corrected 
based on our refined knowledge of the structure forma- 
tion scenarios, based on the perturbation theory as well 
as the N-body simulations |6lL 163 . l65l l6r| . Further- 
more, as an interesting possibility, coupling of different 
Fourier modes induced by non-linear gravitational clus- 
tering could induce a transfer of the dark energy cluster- 
ing effect at large scales to the power of P(k) at smaller 
scales. If this is the case to be apparent on the measured 
power spectrum, an inclusion of the non-linearity in the 
model predictions could improve a signal-to-noise to test 
smoothness of dark energy, from the smaller-scale clus- 
tering information. This issue will be worth exploring in 
detail, and will be presented elsewhere. 
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